PSEUDOCODE 


We present basic pseudo-code for some of the algorithms, 
discussed in the Steps. In our experience, students do benefit 
by studying the pseudo-code of a method at the same time as 
they learn it in a Step. If they are familiar with a programming 
language, they should attempt to convert at least some of the 
pseudo-codes into computer programs, and apply them to the 
set Exercises. 

Bisection Method (STEP 7) 

The equation isf(x) = 0. 


1 read /i, e 

2 repeat 

3 x = (a + t 0/2 

4 if fix) “ 0 then do; 

5 print l Rool is 1 

6 stop 

7 endif 

if /(a) #■ fix) ;> 0 then do: 

9 €4 — X 

1 0 else do: 

J 1 A * j 

12 endif 

1 3 until t — a < e 

1 4 print 'Appro x i m eUion to root ip 1 f je 


Points for study: 

1. What is the input used for? 

2. Explain the purpose of Lines 8 - 12. 

3. Amend the speudo-code, so that the process will always 
stop after preset M iterations. 

4. Amend the pseudo-code so that the process will stop as 
soon as i/wi < s 

5. Write a computer program, based on this speudo-code. 

6. Use your program to solve Exercises 1 and 2 in the Applied 
Exercises. 


Method of False position (STEP 8) 



The equation isf(x) = 0. 


1 lead o, fr, f 

2 reptat 

3 * *= (a * f{b) -b* f{a)yif{b) - /(«)) 

4 if /(*) = 0 then doi 

5 print '“Ronl is', x 

6 stop 

7 endif 

8 if /( B )* /(jt> > Othondni 

9 o = X 

10 else do; 

11 b-x 

12 endif 

J3 until \f(x}\< ( 

14 print ‘Approximation to root is', j r 


Points for study 

1. What are the input values used for? 

2. Under what circumstances may the process stop with a 
large error in x? 

3. Amend the pseudo-code so that the process will stop after 
M iterations, if the condition in Line 13 is not satisfied. 

4. Write a computer Program based on the pseudo-code. 

5. Use your program to solve Exercises 1 and 2 in the Applied 
Exercises . 

• Newton-Raphson iterative method (STEP 10) 

The equation isf(x) = 0. 


1 read lift if 

2 

3 repeal 

4 /> - f (.a)/ f(ti) 

5 a * a — 3 

6 N = M + I 

7 until |*| < t or N = M 

R p lint 'Appro m at i on to root is 1 „ d 
y if Ml > if then do: 

SO prim 'ra|iiiftd accuracy not reached in', , "iterations' 

1 S endif 


Points for study 
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1. How are the input values used? 

2. Why is M given in the output of Line 10? 

3. What happens if f(a) is very small? 

4. Amend the pseudo-code to take suitable action iff (a) is 
very small. 

5. Write a computer program based on the pseudo-code. 

6. Use your program to solve Exercises 1 and 2 in the Applied 
Exercises . 

• Gauss Elimination (STEP 11) 

The system is: 


G|I*I + <712*2+ ■ 

=*l 

021*1 + 022*2 + ‘ 

■ - + = ^2 

0ji!*I + a n2*2 + r 

1 ’ + Jfl — bn 


1 read flu h . , , i fljtRi j!?i , , r , i 

2 for it = 1 ton - 1 do; 

3 fori -k + 1 tond<>: 

4 m = fljjt/djt* 

5 for j = k + 1 to n do: 

6 aij s= Uij - w * fly 

7 endfor 

8 = fij - m * h 

9 endfor 

10 endfor 

11 X n — tft/flnfl 

12 for f = fi — 1 duwnto 1 do: 

13 Xi=ln 

14 for j - i + Hondo: 

15 .Tp - rtf “ HjJ * Xj 

16 tndfor 

17 xi = Jfj/Oif 

IS endfor 

19 print "Approximate solution k\ jti, . - . , 


Points for study 

1. Explain what happens in Lines 2-10. 

2. What process is implemented in Lines 11 - 18 ' ? 

3. Amend the pseudo-code so that the program terminates 
with an informative message when a zero pivot element is 




found. 

4. Write a program based on the pseudo-code. 

5. Use your program to solve Exercises 3 and 4 in the Applied 
Exercises . 

• Gauss-Seidel Iteration CSTEP 13) 

The system is: 


filial + <712*2+ ■ 

- + ai„x„ = b 1 

021*1 + 022*2 + ‘ 

■ ■ + Bln*A a= b 2 

0ji!*I + a n2*2 + ' 

’ ’ + ~ = bn 


1 read fi n aj t fljira 

2 for It = 1 to/i - 1 do; 

3 for i = k + 1 to n d»: 

4 m = flft/djt* 

5 for j = k + 1 to n do: 

6 aij " a/; - w * a $ 

7 endfor 

8 b, = b, - m * bt 

9 endfor 

10 endfor 

11 Xn =^= b^/Onn 

12 for i = n — 1 duvmtu 1 dot 

13 Xi=bi 

14 for jf — i + 1 to/t do: 

15 I, = j, “ H/7 * Xj 

16 endfor 

17 Xi=Xi/Qii 

1 & endfor 

19 print ‘Approximate sotution is*, *i , .... jt« 


Points for study 

1. What is the purpose of the number s? 

2. What are theyi,y2 1> • • • ,yn used for? 

3. Why is it possible to replace theyy in Line 13 by xj? 

4. Amend the pseudo-code to allow a maximum of M 
iterations. 

5. Write a program based on the pseudo-code. 

6. Use the computer program to solve the system: 




8* + y — 2z = 5 

x-7y+z = 9 

2x + 9z= 11 

7. Use your program to solve Exercises 3 and 4 in the Applied 
Exercises . 

• Newton divided difference formula (STEP 24) 

You are to calculate for given data xq, x\, . . . , x n ,f(xo),f(xi ), . . 
.,f(x n ), and given ****>.**]>, the interpolating polynomial P n (x> of 
degree n. (The algorithm is based on divided differences.) 


\ read a, jto t - , . , /(*o) /to* 

2 for [ = 0 to it doi 

3 = /<*> 

4 endfor 

5 for i ■ 1 tos d«t 

G for ; = l to i do: 

8 endfor 

$ endfor 

10 sum^di i.o 

11 prods 1.0 

12 for i m 1 ton do: 

1 3 prod = prod * (je - *<_t } 

14 sum *= jinn ■+■ di, t * prod 

15 endfor 

id print b Appro*inialioji at x *»', jwm 


Points for study 

1. Follow the pseudo-code through with the data n = 2, .x = 
1.5, xo = 0,/Y.xo) = 2.5,.xi = l,f(x±) = 4.7, xSS 2 = 3, and / 
X 2 ) = 3.1. Verify that the values da calculated are the 
divided differences /( xq, . . .,.x/). 

2. What quantity (in algebraic terms) is calculated in Lines 
10-15? 

3. Amend the pseudo-code so that the values Pi( x), 

F> 2 ( x)P(x), • • Pr- i( *)are also printed out. 

4. Write a computer program based on the pseudo-code. 

5. Use your program to estimate f[2) for the data given in 1 
above. 

6. For the data, given in Exercise 6 of the Applied Exercises, 




use the program to obtain an estimate of /o(0.25). 

• Trapezoidal Rule(STEP30) 

The integral is: 




] readd, /?, jV, A/, e 

2 done = false 

3 1/ = 0,0 

4 repeat 

5 ft = (b - a)/N 

6 j = (/(«) +- fib)}/2 

7 for i = 1 to N ~ 1 do: 

8 x = « + i * A 

9 r = j + /(jt) 

1C endfor 

11 r = 

12 if |T - U\ < e then do: 

13 done = true 

14 else do: 

15 N = 2*N 

lfi U = T 

17 end if 

18 until & > M or done 

19 print 'Approximation to integral is*, T 

20 if N > M then do: 

21 print 'required accuracy not reached with M =\ M 

22 endif 


Points for study 

1. What are the input values used for? 

2. What value (in algebraic terms) does Thave after Line 11? 

3. What is the purpose of Lines 12-17? 

4. Write a program based on the pseudo-code. 

5. Apply your program to Exercises 7 and 8 of the Applied 
Exercises . 

• Gauss integration formula (STEP 32) 


The integral is: 






a 


Use the Gauss two-point formula. 


1 rtadd.fr 

2 jej + 

3 XI =» (6 + a + (i - a)}J%)fl 

4 i = (/(*]) + /(«M/a 

5 print 'Approx imaiion to miff ral is’, / 


Points for study 

1. What is the purpose of Lines 2 and 3? 

2. What changes are required to produce an algorithm based 
on the Gauss three-point formula? 

3. Write a computer program based on this pseudo-code. 

4. Use your program to solve Exercises 7 and 8 of the Applied 
Exercises . 

5. Runge-Kutta method (STEP 33) 

Process the equation^’ =f(x,y) and use the usual fourth- 
order method. 


1 nadx,y,fi,M 

2 print x> y 

3 N = 0 

4 repeal 

5 J i y =h*f{x,y) 

6 x^x + h/2 

1 z = y + *i/2 

8 * z = A 

y z -y + fci /2 

10 In a k + /(*. t) 

11 + h/2 

12 + fc} 

13 *4 - h * /(*> *) 

14 y m y + (ftj + + 2*J + 

15 print *. y 

10 .V = N -H 

1 7 until N — M 


Points for study 


1. What are the input values used for? 





2. How many times is the function /evaluated between 
Lines 4 and 17? 

3. Amend the pseudo-code for use with the second-order 
Runge-Kutta method. 

4. Write a computer program based on the pseudo-code. 

5. Use the computer program to solve Exercises 9 and 10 
of the Applied Exercises . 


